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Abstract 

It is known that a semi-classical analysis is not always adequate for atomtronics devices, but that 
a fully quantum analysis is often necessary to make reliable predictions. While small numbers of 
atoms at a small number of sites are tractable using the density matrix, a fully quantum analysis is 
often not straightforward as the system becomes larger. We show that the fully quantum positive-P 
representation is then a viable calculational tool. We postulate an atomtronic phase-gate consisting 
of four wells in a Bose-Hubbard configuration, for which the semi-classical dynamics are controllable 
using the phase of the atomic mode in one of the wells. We show that the quantum predictions of 
the positive-P representation for the performance of this device have little relation to those found 
semi-classically, and that the performance depends markedly on the actual quantum states of the 
initially occupied modes. We find that initial coherent states lead to closest to classical dynamics, 
but that initial Fock states give results that are quite different. A fully quantum analysis also 
opens the door for deeply quantum atomtronics, in which properties such as entanglement and 
EPR (Einstein-Podolsky-Rosen) steering become valuable technical properties of a device. 

PACS numbers: 03.65.Yz, 03.75.Lm, 02.50.Ey, 67.85.Hj 
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I. INTRODUCTION 


Atomtronics is an emerging area of investigation in which analognes of electronic cir- 
cnits and devices are constrncted nsing nltra-cold bosonic atoms rather than electrons as in 


conventional electronics [^. The conventional way to constrnct an atomtronic device is to 
nse cold atoms trapped in an optical lattice, which has a description in terms of either the 
Mott-Hnbbard model transferred from condensed matter physics {2 4| for fermionic atoms, 


fl. 


and the Bose-Hnbbard model for bosonic atoms l5|l. These models can represent either in- 
snlating behavionr, in the Mott insnlator regime, or condncting behavionr, in the snperflnid 
regime. In this work we consider only bosonic atoms. Shortly after the realisation of trapped 
Bose-Einstein condensates (BEC), Jaksch et al. jb] showed that the Bose-Hnbbard model 
can provide an accnrate description of bosonic atoms trapped in a deep optical lattice. The 


basics of this model have been nsed to investigate a wide variety of atomtronic c 


inclnding one with diode-like behavionr j^, a single-atom switching transistor [8|, and one 


proposals inclnde circnits of diodes and transistors 


which nses a triple well conhgnration to mimic a held effect transistor j^. More recent 


B 


evices 


Q. 


lOj, and an atomtronic battery |l1 1. 


Gajdacz et al. proposed atomtronic transistors with the idea of combining them into gates 
for qnantnm compnters 1^. We have no donbt that there will be more proposals in the 
fntnre, and what we will show here is that a fnll qnantnm description of both the dynamics 
and the initial state needs to be taken into acconnt to gnarantee accnrate descriptions of 
the dynamics of snch devices. 

Varions theoretical methods have been used to analyse these devices up until the present. 
For small numbers of atoms and sites, direct quantum calculations using a master equation 
are often possible. Pepino et al. developed a quantum master equation to treat systems 
which interact with sources and sinks, based largely on methods which have been extremely 
successful in quantum optics. With this method, they have analysed electronic diodes, held- 
effect transistors, bipolar junction transistors, and an analogy to a logic gate 13|. The master 
equation was then solved numerically, giving fully quantum solutions. This method allows 
for the investigation of systems with relatively small numbers of atoms and lattice sites, 
with the authors typically treating 3 or 4 sites, each with 1 or 2 atoms. Gajdacz et al. used 
coupled nonlinear Schrodinger equations in one dimension to calculate the phase evolution 
of eigenstates, with two distinguishable atoms in a triple well 12|. For the small systems 


2 












considered by these authors, their methods are perfectly adequate. In the mesoscopic regime, 
when the numbers of atoms can make density matrix calculations complicated, the one- 
body Schrodinger equation can be used for non-interacting atoms, or the mean-held Gross- 


Pitaevskii equations (GPE) for interacting atoms. Bo 


coherent atomic transport in a three-well potential 


15 


h of these have been used to model 


17| and each has disadvantages. The 


one-body Schrodinger equation, being linear, cannot include the interactions required for 
the Mott insulator regime. The GPE approach cannot describe any quantum statistical 
features, such as the actual quantum states or any entanglement properties. The three-well 
coherent transport model has previously been analysed using the fully quantum positive- 
P representation |l8|, which was used to note the differences caused by different initial 
quantum states on the dynamics, the entanglement properties and the quantum steering 
properties [^ . 

In what follows we base ourselves on the approach taken by Milburn et al. 2l|, gen- 
eralisng this to four wells |22|, l23|, and using the fully quantum positive-P phase space 
representation. We consider this to be the most suitable approach here because it is exact, 
allows for an eas y re presentation of mesoscopic numbers of atoms, and can simulate different 


quantum states 


24l |. Just as importantly, the positive-P calculations scale linearly with the 


number of sites and can in principle deal with any number of atoms. These are powerful 
advantages when we wish to consider mesoscopic numbers of atoms in arbitrary numbers 
of potential wells. One disadvantage of the positive-P representation is that the integration 
has a tendency to diverge for high collisional nonlinearities 25|], although it often converges 


for long enough to show marked differences 
for the calculation of quantum correlations 


fom mean-held predictions and also allows 


27|. 


II. PHYSICAL MODEL, HAMILTONIAN AND EQUATIONS OF MOTION 

The four-well system is as shown in the schematic of Fig. [H where the circles represent 
the wells or lattice sites, each of which contains a single atomic mode. The clj are bosonic 
annihilation operators for atoms in mode j, the Ej are the ground state single-atom energies 
of the wells, and J represents the coupling between the wells. We assume that there is no 
coupling between well 4 and wells 1 and 3. We also assume that any atoms at each site 
are initially in their ground state. The basic idea is that atoms from the outside wells will 
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FIG. 1: (Color online) A simplified schematic of the device. Wells 1,2 and 3 have nearest neighbour 
couplings, and well 4 is coupled with well 2. The aj are bosonic annihilation operators for atoms 
in mode j, the Ej are the ground state single-atom energies of the wells, and J represents the 
coupling between the wells. 63 represents a phase-shift which is applied to mode 3 at the beginning 
of the evolution, when the populations of 1 and 3 will be equal, with 2 and 4 initially being empty. 

tunnel into the centre well and from there to well 4. Because of the wave nature of the 
condensed atoms, we expect that the initial phase difference, 63 between the atomic modes 
at sites 1 and 3 will affect the rate of tunnelling into the middle and hence into well 4. We 
would hence have a type of phase sensitive gate, where the phase of one mode mode can be 
used to control the occupation of another. 

Following the usual procedures j^, with y as the s-wave collisional term, we write the 
Hamiltonian as 


n = 



( 1 ) 


where j runs from 1 to 4. Starting from this Hamiltonian, our hrst step is to hnd the semi- 
classical mean-field equations in the Gross-Pitaevskii approach. We will use the solutions of 
these for comparison purposes, since it is well known that they are not always accurate for 
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FIG. 2: (Color online) Solutions of Eq. [2]for the numbers in each well as a function of dimensionless 
time, Jt. The initial phase difference, 63 , between wells 1 and 3 is zero, and the numbers in these 
two wells are equal. Wells 2 and 4 are initially empty. We see that atoms pass through the middle 
well and populate the fourth. The parameters used in this and all subsequent plots are Ej = OVj, 

Ni{0) = NsiO) = 50, A^2(0) = A^4(0) = 0, J = 1, and x = 10“^. 


the Bose-Hubbard model, even for calculation of the mean helds 
aj to represent the mean-helds, we hnd 

dai 


21 


23| . Using the variables 


dt 

d(y.2 

dt 

da^ 

dt 

da^ 

dt 


= -i (Ui + 2x|ai|^) ai + f Ja2, 

= —i {E2 + 2x|tt2p) + id + 03 + 04) , 
= (U3 + 2 x|a 3 |^) 03 + f Jq!2, 

= —i (U4 + 2 x|a 4 p) 04 + f^02- 


( 2 ) 


To solve the full quantum equations, we use the positive-P representation 18|, which 
allows for exact solutions of the dynamics arising from the Hamiltonian of Eq. [1] in terms 
of normally-ordered operator averages, in the limit of the average of an inhnite number 
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of trajectories of stochastic differential equations in a doubled phase-space. This method 


also allows for the representation of different quantum states in the initial conditions 


24|. 


The positive-P representation does have the disadvantage that it can suffer from divergence 
problems in cases with a largish colliisional nonlinearity, often limiting its utility to short 
time dynamics 2^. In this letter, it is exactly this regime we are analysing, so that the 
positive-P representation is perfectly adequate. Following the standard methods Q, the 


set of Ito stochastic differential equations 


28| are found as 


ddi 

dt 

daf 

dt 

da2 

dt 

da2 

dt 

da^ 

dt 

da^ 

dt 

dct^ 

dt 

daf 

dt 


-i [El + 2xafai) ai + iJa2 + \J-2ixoi\ r]i, 
i [El + 2xaiai) af - iJa^ + \j 2ixaf ^ r/ 2 , 

—i [E2 2x0^20:2) C(2 + iJ {oii -|- 0(3 -t- 0:4) -h \J—2ixc^2 V 3 i 

i [E2 + 2xa^a2) - iJ [ctf + \J21x0.2 t] 4 ., 

-i [E3 + 2x0^03) 03 + iJo2 + \J-2ixol r]5, 
i [E3 + ‘ 2 x<^t<^ 3 ) Oit -iJo^ + \J2ixot ^ rjQ, 

-i [Ei + 2xq!4 Q!4) Q!4 -1- iJo 2 + \J- 2 ixol tjy, 
i [E 4 + 2 x 04 04 ) 04 - iJo 2 + \J 2 ixol ^ t] 8 , 


( 3 ) 


where the rjj are standard Gaussian noises with ffj = 0 and 'r}j{t)r}k{t') = 6 jkS[t — t'). As 
always, averages of the positive-P variables represent normally ordered operator moments, 
such that, for example, oj'o'l'^ —)■ 


III. RESULTS 

In order to analyse and compare the effects of different initial states, we begin with a total 
of 100 atoms, evenly distributed between wells 1 and 3 and analyse the classical evolution 
given by Eq. [2J The main quantity of interest is N 4 {t), the number of atoms in well 4 as 
a function of time. With the initial phases equal, we see a 90% transfer of atoms into the 
fourth well at certain times, as shown in Fig. |2l When we set 63 to vr, this fourth well 
remains unpopulated, showing that the device can act semi-classically as a phase-dependent 
gate. With 63 = 7i/2, as seen in Fig. |3l we see that a little less than half the atoms appear 
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FIG. 3: (Color online) Solutions of Eq. [2] for the numbers in each well for the same parameters 
and initial conditions as Fig. [2]except that 6 ^ = 7r/2. We see that the maximum population of well 
four is approximately half that seen in Fig. [2j 


in the fourth well. Having thus benchmarked the device, we will now proceed to see how it 
performs with initial coherent and Fock states in wells 1 and 3. 

The initial quantum states are modelled as in ref. 2^, with the positive-P representation 
equations being numerically integrated over a large number of stochastic trajectories. The 
immediate difference is that, because we can calculate close approximations to normally- 
orderd operator expectation values, we are also able to calculate the standard deviations 
about the average solutions for the atom numbers. We begin with coherent states with the 
same average occupation numbers as used classically. The important differences are that the 
initial physical number distribution is now Poissonian (close to Gaussian as the number of 
quanta increases) and quantum noise is included in the calculations. We immediately hnd 
that the average solutions for the intensities are the same as the semi-classical mean-held 
solutions, for both ^3 = 0 and 9^ = vr. The only difference is that there is a distribution about 
these solutions, as shown by the lines plotted at plus and minus one standard deviation. Thus 
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FIG. 4: (Color online) Quantum solutions for the number of atoms in the fourth well, beginning 
with coherent states, and the mean plus or minus one standard deviation. The initial conditions 
are the same as for Fig. [2] and the solutions are the average of 6.6 x 10^ stochastic trajectories of 
the positive-P equations. 

the semi-classical treatment gives the most likely solutions, although we must remember that 
most individual experiments will give something different. How important this is in the case 
of initial coherent states will depend on how precisely we wish a particular atomtronic device 
to function. 

When we begin the simulations with Fock states of dehnite atom number in wells 1 and 
3, we hnd that the device loses its phase dependence completely. The average results are 
almost identical to those found semi-classically for 6 ^ = 7r/2, and do not change irrespective 
of the initial phase. This is not unexpected since Fock states can be roughly thought of 
as having indeterminate phase. As can be seen in Fig. [6l however, the standard deviation 
around the mean solution is signihcant. In practice this means that the device would be 
totally inappropriate as a phase sensitive gate if the initially occupied wells are in number 
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FIG. 5: (Color online) Quantum solutions for the number of atoms in the first and third wells, 
beginning with coherent states, and the mean plus or minus one standard deviation. and 
are equal. The initial conditions are the same as for Fig. 01 but with 63 = tt, and the solutions are 
the average of 6.45 x 10^ stochastic trajectories of the positive-P equations. 


states. As this is usually thought of as the natural state for Mott insulators, this would be 
problematic. 

Although we have considered the two extreme cases of coherent and Fock states here, it 
is worth noting that there are other possible initial states. As an example, the collisional 
nonlinearity of condensed atoms has been shown to give a “crescent” state, where the Wigner 


function is stretched in the phase space [29j. T 


affect the dynamics of molecular association 


lis initial state has been previously shown to 


30 


31| and would be expected to have an effect 


with our device as well. The advantage of the positive-P representation is that it allows us 
to consider arbitrary initial quantum states, as well as the calculation of any correlations 
that can be expressed as normally-ordered operator expectation values. 
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FIG. 6: (Color online) Quantum solutions for the number of atoms in the fourth well, beginning 
with Fock states, and the mean plus or minus one standard deviation. The initial conditions are 
the same as for Fig. [2] and the solutions are the average of 3.87 x 10^ stochastic trajectories of the 
positive-P equations. The solutions for initial Fock states are identical for all values of 63 . 

IV. CONCLUSIONS AND DISCUSSION 

In conclusion, we have shown that the positive-P representation is a useful tool for 
analysing the performance of atomtronic devices where mesoscopic numbers of sites and 
atoms are involved. It is easily extended to much greater numbers of both than we have 
considered here. We have also shown that mean-held claculations can be very misleading 
for such systems and it is therefore important to undertake full quantum calculations. We 
have considered two markedly different initial states, showing that the coherent states, as 
expected, give results close to the semi-classical predictions. On the other hand, if the atoms 
are initially in Fock states, the semi-classical predictions are almost useless. As shown pre¬ 
viously, correlations between the atomic modes at each site can also be calculated with this 
representation, which will open the door for fully quantum atomtronics in which use can be 
made of properties such as entanglement and EPR steering. This will dehnitely open up 
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regimes that are not readily accessible with standard electronics. 
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